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I INTRODUCTION 

The deep relationship between conformal field theory and criticality has pro- 
vided a wealth of detailed information on phase transitions and critical phe- 
nomena. Moreover, perturbed conformal field theory provides a description of 
the approach to criticality in certain models . One of the most striking ex- 
amples is the (t>(i,2) perturbation of the minimal unitary conformal field theory 
Ms^i which is known to describe the scaling limit of the two-dimensional Ising 
model at r = Tc in a magnetic field. In particular, Zamolodchikov's construc- 
tion of nontrivial local integrals of motion and thus an integrable quantum field 
theory led to the remarkable prediction of eight fundamental mass ratios for 
the magnetic Ising model The masses coincide with the components of the 
Perron- Frobenius vector of the Cartan matrix of the Lie algebra Eg. 

In another development, the exactly solvable dilute A3 lattice model was 
discovered |^ and (in regime 2 of its four regimes) seen to be in the same 
universality class as the magnetic Ising model. Most importantly the dilute 
model ^, ^ admits an off-critical extension in which the Boltzmann weights 
are parametrised in terms of elliptic theta functions In the dilute A3 model 
the elliptic nome plays the role of magnetic field. Its hidden Eg structure has 
been revealed by a number of studies [D-|13 . The masses, obtained from the 
eigenspectrum, may be summarized by the formula ^| 

a \ J / 

where index j labels the eight particles, 5 = 30 is the Coxeter number for Eg 
and the set of allowed a values is given in Table |. 

In addition to the correspondence between the dilute A3 model and Eg, there 
are similar correspondences between the dilute A4 model and E7, and the dilute 
Ag model and Eg. In regime 2 these models are lattice realizations of the 4'{i,2) 
perturbation of the A^4,5 and Ale,? minimal unitary conformal field theories 
respectively, known to have connection to the other exceptional Lie algebras 
[14 . Some E-type structures have been observed for these dilute A models 

Based on the results for the eigenspectrum of the dilute A3 model ||T^ and 
general inversion relations, we proposed [ p7| |l8j that, in the thermodynamic 
limit and in the appropriate regime, the row transfer matrix eigenvalue excita- 
tions 

r,iw) = hm (2) 

of the dilute A3, A4 and Ag models are given by the following general expression. 
Proposition. The excitation spectrum of the dilute A3, and Aq models in 
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regime 2 is given by 



'^jM = ll^ TZ 671^ ■ (3) 

a E{ — X 9 W,X^'^'')E{ — X 9 w,x^^'^) 

Here the elliptic nome is p = e~^, w = 6"^^^"/% and x — e^^^/''^. Regime 2 
is specified by the range of the spectral parameter: < it < 3A, and the value 
of the crossing parameter: A = ns/r where s = L + 2 and r = 4(L + 1). For 
the dilute A4 model the E7 Coxeter number is g — 18, while for the Ag model 
the Eg Coxeter number is g — 12. The standard (conjugate modulus) elliptic 
function is defined by 



E{z,q)^l[{l-q"-'z){l~q"/z)il 



n=l 



The numbers a appearing in (||) are given in Tables |, || and III . The integers 
in these tables have appeared in other contexts in relation to the E-algebras 

In this paper we explicitly derive the elementary excitation spectrum of the 
dilute A4 model, thereby confirming our Proposition in this case. The result 
(||) leads to the inverse correlation lengths and mass gaps. Our input to these 
calculations are the string solutions to the Bethe equations found by Grimm 



and Nienhuis ||9|, |10|, |21|]. As discussed later in IV, our results are applicable to 
the tricritical Ising model which is in the same universality class. In particular, 
the elliptic nome appearing in the dilute A4 weights in regime 2 corresponds 
to the leading thermal off-critical perturbation in the tricritical Ising model. 
This perturbation is identified with (/)(i,2) and has been shown to exhibit 
E7 structures We are able to obtain exact results for some universal 

amplitudes of the tricritical Ising model. These results are in agreement with 
those found recently by other means ^ . 

The outline of the paper is as follows. The dilute Al lattice model is defined 
along with the corresponding Bethe equations in The bulk free energy and 
the eigenvalue expressions in regime 2 for L — A associated with the seven E7 



masses are derived via the exact perturbation approach in III (continued in the 
Appendix). The paper concludes in ^ with a discussion of the results and their 
relevance to universal behaviour in the tricritical Ising model. 



II THE DILUTE A. MODEL 



We here give a short summary of facts about the dilute A^ models 26, g_3| 
which are pertinent to our calculations. 

The dilute Al model is an exactly solvable, L-state restricted solid-on-solid 
model defined on the square lattice. Its adjacency diagram is the Dynkin dia- 
gram of Al with the additional possibility that a state may be adjacent to itself 
on the lattice. The model is solvable in four off-critical regimes, with the elliptic 
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nome p of its Boltzmann weights taking the model off-critical. At criticality, 
the dilute model can be constructed Q from the dilute 0(n) loop model 
[p7| p8[ . In regime 2 of the model the central charge is 
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L{L + iy 



In the majority of exactly solved models the elliptic nome plays the role of 
temperature |p9| . In the dilute A^ model the interpretation of the elliptic nome 
differs according to whether L is even or odd. For L odd the elliptic nome plays 
the role of a magnetic field [||, and p > and p < are related by simple 
label reversal of the heights. For L even the nome plays a thermal role, and the 
behaviour of the model depends on whether p > Q (regime 2"*") or p < (regime 
2~). More specifically, it was shown |^ that in regime 2 the nome corresponds 
to perturbation of the A^l.l+i minimal unitary conformal field theories by the 
operator c^(i,2)- 

Using the conjugate variables introduced after (^), and setting Wj — e^^Truj/e^ 
the eigenvalues of the row transfer matrix of the dilute A models (for a peri- 
odic strip of width N where for convenience N has been taken as even) can be 
written js) 



k{w) =UJ 



£:(a;4«,a:2'-) E{x^'' ,x'^'^) 
a;^'* E{w,x'^'') E{x^' /w, x^'') 



N N 



n 

N 



l-lsjr E{x^''wlWj,X^'') 

^ Eix^^Wj /w,x'^'^) 



w i;(a;4^a;2'-) £'(x6^a;2'') 

N 

X w 
J 

„2r\ T?( ^2s /„,, _2r 



E{w/wj,x'^'') E(x^^Wj/w,x'^'') 
\ "'^ E{x'^''Wj/w,x'^'') E{x'^''Wj/w,x'^'') 



+ LU~ 



2^ E{w,x^'') E{x^' /w,x^'') 
E{x^^,x^-') E{x^^,x^'-) 



N N 

n 



w 



2s/rE{x^''Wj/w,X^'') 



E{x'^^Wj /w, X 



2r^ 



(4) 



where uj = exp{in£/{L -I- 1)) for ^ = 1, ... , L, and s = L + 2 and r = 4(L + 1) 
in regime 2. The Bethe equations which give the N roots Uj have the form 



Eix^^/wj^x^"-) 



-, N 



E{x 



2s, 



„2r 



N 

n 

k=l 



2s/r Ejx'^^Wj/wk, X^"") E{x'^''Wk/Wj,x'^'') 

'^^ E{x'^^Wk/wj,x^'-)E{x'^'^Wj/wk,x'^'-)' ^ ' 



In the limit \p\ — > 1 with u/e fixed, or equivalently x —f 0, the excitations in 
the eigenspectrum Vj (w), defined in (^), break up into a number of distinct bands 
labelled by integer powers of w. Numerical investigations of the eigenspectrum 
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[|[ ^, |l^, ^ have revealed eight and seven thermodynamically significant 
excitations for L = 3 and L = A respectively, and provided the data in Table 

m 

We previously |l^ applied the exact perturbation approach initiated by 
Baxter |Q to calculate the excitations in the eigenspectrum for L = 3. This 
involved perturbing away from the strong magnetic field limit at p — s- 1; for 
L = 4 this limit corresponds to moving far away from the critical temperature. 
The calculations follow. 



Ill MASS SPECTRUM 

A Preliminaries 

To apply the perturbation technique pO{ to find the form of the excitations 
the string structure of the Bethe ansatz roots (||) is required input. The 
groundstate roots all have uj pure imaginary, so that for 
j = 1^... ,N with \aj\ = 1; in this sense they all live on a unit circle. For 
each excitation i, certain roots acquire a real part m7r/20, as shown in Table 



[V . (If there are rii such roots, one says there is an rii-string associated with 
the excitation.) For these roots Wj = hjX^ , so that the string entries can be 
thought of as living on circles of radius a;"* with phase hj , while the other N — rii 
roots again lie on the unit circle. 

The process of finding the excitations involves using the Bethe equations 
to set up recurrence relations for auxiliary functions of the unknown roots 
Qj. As the roots only enter the eigenvalue expression (^) through the auxiliary 
functions, it just remains to solve the recurrence relations by iteration and to 
simplify the resulting expressions. The largest eigenvalue Ag, relative to which 
excitations are measured, was calculated previously in this way for all L. 

The relationship between the excitations (Q), the correlation lengths and 
the mass spectrum nij of the associated field theory is 

= -logrj ^iTij, (6) 

where we take the isotropic value u — 3A/2. 

It is convenient to use the notation for products: 

oo 

{z;pi,... ,pk)oo= Y[ 'y^^PT ---pTz) 

ni,...nfc=0 
m 

(zi, . . . . . . ,pfe)oo = nC^^jSPi, ■ ■ • ,Pfe)oo 
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which satisfy many identities, the ones used repeatedly in what follows being: 

7 ^ = (2:;g)oo 

{zq/p;p, q)oo ^ {zq/p;q)^ 

{z;p,q)oo {z;p)oc 

The standard elliptic function is thus rc-written as 

oo 

E{z, q)^l[{l- - q^/z){l - g") = (z, q/z, q; q)^. (7) 

ri=l 

It also proves convenient to use the shorthand notation IlJLi '^j = 

For each rrii, if the associated string of excited roots has length n^, we define 
the required auxiliary functions of the as-yet-unknown roots to be 

Fi{w) = ]^ (w/aj;a;^'')oo, 

Gril/w)^ n (x2-a,/iz;;x2-)oo. (8) 

In fact, we actually solve for combinations of these: 

= F,{w)/F,{x'''w) = F,{w)/F,{x^--^'w), 
g,{l/w) = G,{l/w)/G,{l/x'^w) = ai/wj/ai/x^^-^'w), (9) 

for i — 2,A, 5, 6, 7 (but for i = 1, 3 slightly different definitions are convenient 
and are given as required). 

So far as possible, we write factors and powers which are common to all 
eigenvalues (or indeed to all the eigenvalues for other models) in terms 
of the generic r and s to distinguish them from the particular integers which 
arise from the input strings. Of course, r — 20 and s — 6 throughout. Once 
the particular string form for the roots has been applied, the calculations are 
straightforward for all masses except mi and 1113. For this reason, we sketch 
below the details for the first three masses. The other cases follow similar paths 
to 7712 or indeed to most of the masses for the dilute A3 model , so we relegate 
them to the appendix. We make some comments concerning mi, and mg 
later on. 

B Mass mi 

We begin the perturbation argument with the structure Wj — aj for j = 
1, . . . , iV — 3 with wn-1 = bix~'^, = 62^;^ and wn — b^x^^, so that 
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the string length is ni ~ 3. From the Bethe equations (||) for j = N — 2, 
j = N — 1 and j = in the hmit x ^ we can show that 6i = ^2 = ^3 = b. 
The Bethe equation for the other roots = a is then 



N 



E{x%/a)E{x^%/a)E{x''^h/a) ^yr E{x^' a/ aj)E{x^' a, / a) 
^ E{x^a/b)E{x'^^a/h)E{x'^^a/h) j}^ E{x'^''aj /a)E{x^''a/aj)' 

In the X — > hmit this gives the equation 



(10) 



(11) 



which is an equation of order {N — 2), so that there is a missing root on the unit 
circle, a 'hole', which we call aN-2- Since this is an equation for the roots, its 
left hand side must be equivalent to IljLi^C'^ ~ ^j): ^^'^ equating the constant 
terms from these two expressions we obtain 



-(AAr_3&3)3/5 ^ AN^2b^ = AN^3aN-2b^, 



(12) 



(which we later apply to prefactors in Ai). The Bethe equations for b taken 
together in this limit, and combined with ( p^ ) give 



We use this, together with the fact that each root Uj, including the hole, must 
satisfy (O), to show 



{aN~2)'^ ^ = - A]^s ^za]^s ^2 =^ (aAf~2)^ = 1- 
We define the following auxiliary functions of the roots (see (^): 

Fiiw) ix^w/b;x^noo 



(13) 



:Fi{w) 
Giil/w) 



Fi{x^^w) {x^^w/b; x^^)oo' 
Gi{l/w) {x^%/w;x^'-)c 



Gi(l/xi6w;) {x^%/w;X^^)oo' 
They must satisfy recurrence relations arising from (^0|) 

ix^^a;x^no 



^i(a) 
Giil/a) 



^ {x^^a/aN-2, a;28a/ajv-2; x^^oo T^{x^'a) 



(2;2r 2sQ,.2;2r~)^ 
(^2r+2./^.^2r)^ 



{x^'^a/aN-2,x'^^a/aN-2] x^'')oc ^i(a;''*a) ' 



{x^''/a;x^'-)c 



{x'^aN-2la,x^^aN-2la\x'^)oogi{x''la) 
{x'^^aN-2/a, x^'^aN-2/a; x2'')oo Gi (x^^/a) 



(14) 
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Solving these we obtain 



ei(l/a) = ao(l/a) 



{x'^^a/aN^2;x'^'')oo {x^^a/aN^2,x'^^a/aN-2; x''-^^)co ' 
{x'^°aN-2/a;x'^'')oo {x^^aN-2/a,x^^aN-2/a;x^'^'')c 



{x^^aN-2/a;x^'')oo ix^°aN-2/a, x'^^aN~2/a; x^^'')^ 



(15) 



Here J-q and Qo arise from the square bracketed factors in ([l^ ) and give rise 
to the square bracketed factor in (^. They are related to the groundstate 
eigenvalue Aq, they are common to the calculation of each mass and we will 
suppress these factors for m2, . . . , m.7. We now write the eigenvalue expression 
in terms of the auxiliary functions, the first term being 



Ai w 

3 aN^2 



^^2r 6s^^^2r ^^yj ^ x^"^ / W , X^^ / W] x'^'') c 



{X^^ W I ttM -2, X^'^ aN ~2 / W] X^'')e 

{x^'^w/aN-2-, x^^aN-2/w;x'^'')c 



'-Ti{x^'w)gi{l/x'''w). (16) 



Substituting the solutions (^_5|) gives an expression for the excitation ri (w) which 
may be written in elliptic functions (^ as 



Ai 
Ao 



w 



E{-x^^/w, x^^')E{-x^^w, x^^') 
Ei^x^'^w, x^'^'')E{-x^^ /w, a;i2«) ' 



(17) 



where we have set aN-2 = — 1- (The other two terms in the eigenvalue always 
give identical elliptic function expressions to the first, upon simplification.) 

The Bethe equations involving b and the 'hole' equation, which is (|l0| ) with 
a — aAr_2, can also be expressed in terms of the auxiliary functions. Application 
of identities and simplification gives: 

E{x^H/aN-2,x^'-^') = E{x^^aN-2/b,x^'^-^') 



E{x^^aN-2,x^^-)E{x^yaN-2,x^^^) 



E{x^yaN-2, a;i2«)£;(a;48a^_2, a;i2^) 



(oat- 



Clearly = 6 = — 1 (identified initially from numerical studies) satisfy 

these conditions; the second reduces to ( p^ ) in the x —>■ limit, and note the 
similarities with (Ki 



C Mass m2 

We begin the perturbation argument with the structure Wj = aj for j = 
1, ... ,7V — 2 with wn-i = bix~^'^ and wn = 622;^^, so that n2 — 2. From 
the Bethe equations for j = N — 1 and j = N we can show that 61 = 62 = b. 



8 



The Bethe equation for the other roots = a is then 



N 



E{x^%/a)E{x^%/a) E{x^'a/aj)E{x^'aj /a) 



(18) 



E{x^Oa/b)E{x^'ia/b) JLJL E{x^'a.i /a)E{x^'>a/aj)' 
In the a; ^ hmit this gives the equation 

which has the same order as the number of unknown roots (N — 2) so that there 
is no hole. Equating this with n^L^^('* ~ ^j) '^^ obtain 



N~2 



(which we later apply to prefactors in A2). From the other Bethe equations in 
this limit, 



-{An -2 



^2)3/5 



{AN-2b'f 



b^^ = 1. 



(19) 



Treating the Bethe equation ( |18D as before gives, in terms of the functions 
defined in (||) and (^), the recurrences 



^2 (a) 
^2(1/0) = 



{x^^a/b, x^'^a/b; x^Qop T2{x^''a) 
{x^^a/b, x^^a/b; x'^'')^ Ti^x^'a) ' 
{x''%la,x''%la-x^^)^g2{x'' la) 



{x^^bja, x^^la; x2'-)^ Q2{x^'la) ' 
Solving these we obtain 

(x^Ofl/fe, x^'^ajb-, x^'-)^ (a;26a/6, x^^a/b, x^^ajb, x^'^ajb; x^'^')^ 



^2 (a) 
^2(1/0) 



(xi4a/6, x26a/5; x2'-)oo (xi^a/b, x'^'^ajb, x^'^a/b, x'^^a/b; x^^')c 
ix^%/a, x''%/a; x^Qop {x^^b/a, x^%/a, x^%la, x'^^/a; x^^-')c 
(x54&/a, x^%/a; x'^^)oo (x^%/a, x^^b/a, x™6/a, x^^/a; x^^^)r_ 



We now substitute these into the eigenvalue expression, the first term of which 
is 

A2 {x'^^w/b,x^^w/b,x^b/w,x^^b/w;x'^^')c 



6^ {x'^w/b, x^^w/b, x'^^b/w, x^^b/w; x 



2r\ 



-J^2{x''w)g2il/x''w). 



This gives an expression for the excitation in elliptic functions (setting b = — 1): 



A2 _ 2 E{-x^/w, x^^') E{-x^'^/w, x^^'')E{-x^^ w, x^^") E{-x^°w, x^^') 
A^~^ E{-x^ w, xi2^) E{~x^^ w, xi2^)i;(-x38 /w, xi2«) E{-x^" /w, x^^^) ■ 



(20) 
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If the product of the six Bethe equations involving b is expressed in terms 
of the auxihary functions, the equation for b (generaHzing = 1 seen in the 
a: ^ hmit in (^9|)) is clearly satisfied by 6 = — 1: 

' E{x^b,x^^') E{x^%x^^')E{x'^yb,x^^') E{x^° /b,x^^') '\ ^ _ 
_E{x^/b, xi2'*) E{x^^/b, a;i2«)i;(x386, a;i2«) E{x^°b, x^^'')\ 

Compare the pattern of powers of x in this equation with those in (|20|); this 
equation has a precise analogue for each mass , . . . , mi , which will not be 
given. 



D Mass rris 

We begin the perturbation argument with the string structure Wj — aj for 



bix ^"^ ^ wn^i ~ b2X^^ and wn 



bsx^". From 



J = 1, . . . , — 3 with 
the Bethe equations for j = N —2 and j = N —1 we can show that 6i = 62 = ct, 
but the Bethe equation for j = N does not link 63 = & to a in the x — > limit. 
(This feature was observed also in the L = 3 case, for a string of odd length 
[O.) The Bethe equation for the other roots = a is then 



E{x^'/a) 



E{x'^''a) 



E{x%/a)E{x%/a) 
7^ E{x^a/b)E{x^a/b) 



E{x^^a/a)E{x^^a/a)E{x^^a/a) W E{x^'a/aj)E{x'^''aj /a) 
E{x^^a/a)E{x^^a/a)E{x^^a/a) E{x^'aj /a)E{x^'>a/aj) 



(21) 



In the a; — > limit this gives the equation 



1 



- -{AN-sa'by^yab' ^0. 

N-3, 



Equating this as usual with Jli=i ~ ^j); obtain 



AN^^ab 



(which we later apply to prefactors in A3). From the other Bethe equations in 
this limit, 



-{AN-3a- 



'■bf"> 



2\3 



(AN-zah^) 



= 1. 



In this case it is convenient to define 

F^{w) (x^'^wja-.x^^') 



F3{x^^w) [x^wja^x^-') ' 

G-i(\lw) (a:28a/w;a;2'-) 
G3(l/2;i6w) (a;36a/w;x2'-)' 



(22) 
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because this choice wih make it clear that a is a spectator in the solution to the 
recurrence relation; it does not appear in the eigenvalue expression. 
Treating the Bethe equation (|2^) as before gives the recurrences 



^^"^ {x^a/b,x^a/b-x^^)oo T^ix^'aY 
, I , (x566/a,a;6"6/a;x2'-)^g3(x^Vfl) 

Solving these we obtain 

(a;36a/6; x''")^ [x'^'^ajb, x^'^a/b, x^^a/b, x^^a/b; x^^')c 



Ma) 
Gsil/a) = 



{x^Oa/b;x^^)oo {x^a/b,x^a/b,x^<ia/b,x^^a/b;x^^^)oc ' 
ix^%/a; x^'')oo (x'Sfe/g, x^^b/a, x^%/a, x^^fc/g; x^^')oo 
{x^%/a; x^^)oo {x^^b/a, x^^b/a, x^»b/a, x'^%/a; ' 



We now substitute these into the eigenvalue expression, the first term of which 
is, in terms of the functions (|2^), 

A3 {x^'^w/b, x%lw;x'^-')o^ 2s (M 2s ^ 

Y = V {x^w/b,x^%/w;x^no. ^'^" ""^^'^ '^ 
With 6 = — 1 this gives the expression in elliptic functions 

A3 _ 2 E{-x^ /w, x^^') E{-x^^ /w, x^^')E{~x'^^ w, x^^') E{-x^^^ , w, x^^') 
A^'"^ E{-x^ w, ^12'') Ei-x^'^ w, 2;i2^)i;(-a;44 /w, x^^'^) E{-x^^ /w, x^^') ' 

(23) 

The Bethe equations involving a and b, also expressed in terms of the auxiliary 
functions, give: 



E{x^H/a,x'''-^') = E{x^^a/b,x''"^') 



E{x%, x^^'')E{x^H, x^'^')E{x^^lb, x^'^'')E{x^^/b, x^^^'Y ^ 



E{x» /b, x^^^)E{x^^/b, x^'^')E{x^^b, x^^'')E{x5'^b, a;i2-) 

Notice that b = —1 satisfies this second equation, and the (so far missing) link 
between a and b is provided by the first. 

E Comments on the 'odd' strings 

To close this rather technical section of the paper, we wish to briefly comment 



on the strings of odd length (see Table IV ) 



For this model strings of odd length appear for the first, third and sixth 
masses. In the first case, the odd string of excited roots is accompanied by a 
'hole' among the roots on the unit circle; it is only this hole an -2 which appears 
in the eigenvalue expression. In the case of the third mass, the phase b of the 
string entry m = 10 appears alone in the eigenvalue expression, and the other 



11 



entries of the string seem to have a spectator role. (For the sixth mass, there 
was nothing special about the calculation.) In the calculations for dilute A3 
[ p^ odd strings were involved for masses 4 and 6, where again the calculation 
of the associated excitation was less straightforward than for even strings. In 
one case, both the coefficient of the 'odd' entry and a hole appeared in the 
eigenvalue, while the other calculation resembles that of TO3 in this paper. In 
general the string entries come in pairs ±m, except for m = r/2, which stands 
alone if it occurs, due to the period of the original elliptic functions in Uj . This 
is the only source of strings of odd length; we can only conclude that when such 
an entry occurs, it in some sense dominates calculations following the exact 
perturbation technique. For strings of even length, all the excited roots seem 
to contribute in a more equal fashion to the calculation and to the resulting 
eigenvalue expression. 



IV DISCUSSION 

In this paper we have made use of the Bethe Ansatz string solutions found by 
Grimm and Nienhuis to derive the excitation spectrum of the dilute A4 model 
via the exact perturbation approach. Our expressions for the seven thermody- 
namically significant excitations for the dilute A4 eigenspectrum in regime 2~ 



are given in (gOD, (|2^), (|Ag), (|AJ), (|A6D and AAT^. In this way we have 



verified for a second case the Proposition given by (H). 

It is perhaps unsatisfying that an elegant closed form expression such as (||) 
has been confirmed in the A4 case by relying on numerical data for the strings 



(Table IV). Indeed, as described for the A3 case in |13 , and in the detailed 
study |lC| , tracing the strings from p = Q (criticality) to the position they take 
in the scaling (massive) limit reveals complicated structure (reported with one 
difference by two groups of authors 0). Fortunately, (||) was conjectured 
[ p7| on the basis of general properties of the dilute A models and of the 
E-type algebras, known to be linked by their common connection to the ^(1,2) 
perturbation of the minimal unitary series; the (scaling limit) string data used 
here has not contradicted it, and (admittedly limited) numerical studies agreed 
with the lower eigenvalues jl^. A forthcoming paper Q should shed some 
new light, from the perspective of Coxeter geometry, on the excitations (|^) and 
hence, among other things, on the string conjectures to which they are related 
as demonstrated here in the L = 4 case. 

Recall that the central charge for dilute A4 in regime 2 is c = jq- There 
are several other known manifestations of the c = theory. The Blume-Capel 
model is related to the Blume-Emery-Griffiths model |Q, a classical spin-1 
Ising model with lattice Hamiltonian 

'^BEG = — J ^ SiSj — 15^(1 — Sf) — h''^ Si — H3 ^ SiSj{Si + Sj), 

(24) 

where J is the nearest-neighbour interaction, £> is a crystal field, H a magnetic 
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field term and H3 is a staggered magnetic field. The phase diagrams of these 
models exhibit a tricritical point, as had been observed in physical systems |Q. 

The critical exponents, known from renormalization group studies, are re- 
lated to the Kac table of the c — conformal field theory Q . 

After the Ising critical point, the universality class of the tricritical Ising 
model corresponds to the second simplest unitary conformal field theory in two 
dimensions. It is also the first of the super-conformal minimal models. It can 
be perturbed by its four relevant scaling fields, shown in Table ordered ac- 
cording to the associated conformal weight. The leading magnetic perturbation 
is believed to be non-integrable |Q , and each of the other three perturbations 
give integrable quantum field theories. In the scaling limit these can each be 
associated with a solvable interaction round a face (IRF) model (or to the terms 
in (p^)). The ABF A4 model in regime III ^ realizes the subleading ther- 
mal perturbation. A lattice realization of the subleading magnetic perturbation 
is given by the dilute A3 model in regime 1 and the scaling limit of the 
leading thermal perturbation corresponds to the dilute A4 model as considered 
in this paper. 

The leading thermal perturbation is known to be integrable and massive, the 



masses being described by E7 Toda field theory ]14 



p3| . Numerical results from 



a finite-size analysis in the spin-chain formulation [ 38| , and from field theory via 
the truncated conformal space approach |^ demonstrated the first few masses. 
These are: 

mi = 1 odd 



m2 = 2cosf| = 1.285 575. 

m3 = 2cosf = 1.879 385. 

m4 = 2cos^ = 1.969 615. 

mg = 4 cos ^ cos fl = 2.532 088. 

mg = 4 cos f cos ^ = 2.879 385 . 

m7 = 4 cos cos I = 3.701 666. 



even 
odd 

even (25) 

even 

odd 

even 



The mass spectrum can be classified |g2[ into even and odd states (as in- 
dicated in (p5|)) corresponding to the Z2 symmetry of the affine E7 Dynkin 
diagram. Each of the above seven masses appears in the high-temperature 
phase of the tricritical Ising model. However, only the even subset appears in 
the low-temperature phase. This is consistent with the numerical observations 
on the eigenspectrum of the dilute A4 model 0. For regime 2+, in a study 
of the low- lying excitations, the first and third were absent. As we have demon- 
strated, all seven excitations are present in regime 2~ which (through a quirk 
in labeUing) corresponds to T > Tc- 

Our expression (|^) gives the correlation lengths and related masses (|^), ex- 
pressed in terms of standard elliptic functions and the original nome p, as 

U4(,3g 4,P ' j 



13 



In the critical limit p — > the leading order behaviour is 



5p5/9^sinff. 



(27) 



Substituting the integers of Table applying trigonometric identities and tak- 
ing mass ratios it was demonstrated that the E7 mass spectrum (|2^) is 
recovered. 

The ground states of the tricritical Ising model (in zero magnetic field) have 
been identified |3|, For T < Tc, the system is in a two-phase region of spon- 
taneously broken spin reversal symmetry, with two degenerate ground states in 
the thermodynamic limit. For T > Tc there is one ground state. This ground 
state picture is also consistent with that of the dilute A4 model as \p\ 1. 
In regime 2+ there are two possible ferromagnetic ground states, while in regime 
2~ there is a single disordered ground state. (It is the presence of such disor- 
dered states for L even which complicates the calculation of order parameters 
for this half of the dilute hierarchy.) 

Very recently, an array of universal ratios for the critical amplitudes of the 
tricritical Ising model have been calculated ||2j, |2^ by field theoretic methods. 
Not all of these quantities appear to be accessible via the dilute A4 model. 
However, one such ratio involves the correlation length prefactors ^q, above 
and below the critical temperature. Our results and observations on the eigen- 
spectrum of dilute A4 give this same value: 



— — — — ^ 

We previously jl^ derived the amplitude 

1 



2 cos 



57r 



8%/3cos(27r/9) 



= 0.09420 . 



(28) 



where /s is the singular part of the free energy. This agrees with the determina- 
tion of this quantity for the 4'{i,2) perturbation of the c = field theory Q. A 
related universal quantity is the amplitude ratio associated with the correlation 
length |l| 

where A/ a is the amplitude of the specific heat and a is the related critical 
exponent. Our expressions for these quantities are 
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93V3cos(27r/9) 



0.101678. 



239273 cos(57r/ 18) sin(57r/9) 



= 0.083889 . 
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which agree with the numerical values of p5[ (allowing for a difference in defi- 
nition by a factor a^l"^ ). As remarked p^, ^5|, such values may be observed in 
experimental systems within the tricritical Ising universality class. 
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Appendix A Further mass calculations 

1 Mass m4 

We begin the perturbation argument with the structure Wj — aj for j = 
1, . . . , iV — 4 with wn-3 = bix^^^, — b2X^^ , wn~i = b^x~^ and wn = 

hiX^ . From the Bethe equations for j = — 3, . . . , TV we can show that 
hi — 1)2 — bj, = hi = h. The Bethe equation for the other roots is 



- UJ 



E{x^'/a) 



Eix^'a) 



N 



, 3^sa4 E{x%la)E{x%la) E{x^^ a/ a^)E{x^^ a, / a) 
y ^-4 ) E{x^a/b)E{x<^a/b) jj-^ E{x^''aj/a)E{x'^'>a/aj)' 

In the a; — ^ hmit this gives the equation 

UJ 

so that as usual we find expression involving the prefactors 



(Al) 



UJ 



Using this with the other Bethe equations in the a; — > limit we obtain 



From (Al) come the recurrences 

(x34a/6, x'^^a/b- x^"-)^ T^x^'a) 



Gi{l/a) 

The solutions are 



Qi{lla) 



{x^a)/b,x^a/b;x^'')oo T^ix^'a)' 
{x^%/a, x^%/a; x^^)^ Giix^'/a) 
{x^%/a, x^%/a; x^'-)^^ G^x'^'/a) ' 

{x^>^a/b, x^^a/b, x^°a/b, x^^a/b; x'^'^)^ 
{x^a/b, x^a/b, x^^a/b, x^^a/b; x^"-)^ 

{x^^a/b, x^^a/b, x^'^a/b, x^°a/b\ x^"^')^ 
^ {x™a/b, x'^^a/b, x^^a/b, x^<^a/b; x^^')oo ' 
ix^%/a, x^°b/a, x^%/a, x^%/a; x^^)^^ 
^ (x626/a, x^%/a, x'^%/a, x'^^b/a; x^'-)oo 
{x^%/a, x'^^/a, x^°%/a, x""6/a; x^^')r_ 
^ {x^H/a, x^^b/a, x™b/a, x'^%/a; xi2^)oo 
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In terms of these functions the eigenvalue may be represented as 

3 62 (a;i0u;/6^^22y,/5^2;306/«^,xi85/u;;:E2r)^-^4ia; wjy^y^ fx w). 
Thus, appHcation of the perturbation argument yields the excitation to be 



£;(-.t17w;, E{-x'^lw, x^^')E{-x'^'^w, x^^') E{-x^°w, x^^') 



A4^ 2 

Ao E{-x^°w,x^^'')E{~x^^w,x^^'')E{-x^<^/w.x^^^)E{-x5°/w,x^^'')' 

(A2) 



where we have put b — —1. 



2 Mass ms 



We begin the perturbation argument with Wj = aj for j ~ 1, . . . . iV — 4 and 



Wi\r-3 = &ia;~^^, 'Wn^2 — b2x'^'^, wjv-i — b^x^'^'^, wpf — b^x^'^. We can show 
that the 6,; are equal, and we call them b. The Bethe equation for the other 
roots is 



E{x^''/a) 



E{x'^''a) 



E{x%/a)E^{x^H/a) 
¥ E{x^a/b)E^{x^'^a/b) 



E{x^%/a) E{x^'a/a,)E{x^'a,/a) 



E{x^^a/b) jj-^ E{x^'aj/a)E{x'^''a/aj) 
In the a; — * limit this gives the equation 

id 

which leads in the usual way to a prefactor expression 



From this and the other Bethe equations 



-{A^-^by/^ 



l,5N 



,5N 



1. 



Rearranging (A3), the auxiliary functions obey the recurrences 
(x24a/6, x^Vfc, x^^a/b, x^^a/b; x^'')oo T^{x^'a) 



^5 (a) 



[x^a/b, x^^a/b, x^^a/b, x^^a/b; x^'-)oo J'^ix^'a) ' 



(A3) 



^5(1/0) 



jx^^b/a, x^%/a, x^%/a, x^%/a; x^'')^ Gsix'^a) 
(x486/a, x^%/a, x^^/a, x^%/a; x^^)oo G5{x^''/a) ' 
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The solutions are 



^5(1/0) 



(a;i2a/6, x^^a/h, x'^^ajh; x^^)oo 

jx^'^a/b, x^'^a/b, x^°a/b, x^^a/b, x^'^a/b, x^'^a/b; x^^')^ 
^ {x»a/b,x^'^a/b,x'^Oa/b,x^^a/b,x'^'^a/b,x'^^a/b;x^'^^)^ 
{x^%/a,x'^%/a,x^%/a-x^'^)^ 
' {x^%/a, x'^^b/a, x^%/a; x'^'-)oo 

{x^H/a, x^%/a, x^%/a, x'^%/a, x''%/a, x^'^%/a; x^^')o, 
^ {x^H/a, x^%/a, x^%/a, x^%/a, x'^H/a, x^%/a; x^^')oo 

which we next substitute into the eigenvalue expression 

As _ {x^^w/b, x^^w/b, x^^w/b, x%/w, x^%/w, x^%/w; x^Qoq 
3 b^ {x^'w/b,x^'^w/b,x^^w/b,x'^'^b/w,x^^b/w,x^^b/w;x'^^)oo 

X J^5{x'''w)g5{l/X^'w), 

to obtain (with b = —1) an expression in elliptic functions of nome x^'^^ 



= w 



^E{-x^ /w)E{-x^^ /w)E{-x'^^ /w)E{-x^^w)E{-x'^^w)E{-x^^w) 
E{-x^w)E{-x^'^w)E{-x^^w)E{-x^^ /w)E{-x'^^ /w)E{-x^'^ /w) 



(A4) 



3 Mass me 



We begin the perturbation argument with Wj = aj for j = 1, . . . ,N — 5 and 
wn-4 = bix'^^, wn-3 = b2X~^^, wn-2 = b^x^^ , wn-1 = bix~^, wn = 652;^. We 
can show that the 6, are equal, and we call them b. The Bethe equation for the 
other roots is 



E{x'^'/a) 



N 



{A 



N-5 



^5,3;, E{x^/a)E{x^/a) 



65 E{x^a/b)E^{x»a/b) 

E{x^^b/a)E{x^%/a) E{x^'a/a,)E{x^'aj/a) 
E{x^'^a/b)E{x^^a/b) E{x'^^aj/a)E{x^''a/aj)' 



(A5) 



In the a; ^ limit this gives the equation 



which leads in the usual way to the expression 



-{An -5 



b'f/' = AN_,b'. 



From this and the other Bethe equations 

(A ft5x3/5]'_ {AN-5b^f 

JAN-,b ) - —^w— 



h5N 



1. 
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After rearranging ( [A5| ), the auxihary functions obey the recurrences 
{x^^a/h, x'^^a/h, a;32a/6, ^36^/^. ^2^-,^ Tf^ix'^'a) 



The solutions are 



{x^a/b, x»a/b, x^'^a/b, x^^a/b; a:2'-)^ Teix^'^a) ' 
{x^%/a, x'^%/a, x^%/a, x'^%/a; x^Qoo Geix^Va) 
(x48&/a, x^^b/a, x^%/a, x^%/a; x2'-)^ g&{x^^la) ' 



{x-'^a/b,x^^a/b;x''-)^ 
' {x^<^a/b,x^°a/b; x^'-)oo 

(x^^g/b, x36a/6, x^°a/b, x^^a/b, x^^a/b, x^'^ajb; x^^'), 
^ (x-^ajh, x»a/b, x^'^a/b, x^^a/b, x^^^a/b, x'^'^b/a; a;i2'^)oc 
{x^''b/a,x^%/a-x^'-)oo 
''{x^'%/a,x<^^b/a;x^'-)oo 

{x'^^b/a, x3^5/a, x'^%/a, x^%/a, x'^^/a, a^^^b/g; x^^'), 
(x'>%/a, x^%/a, x^%la, x^Hja, x'^%/a, x'^%/a; x'^^'), 

which we next substitute into the eigenvalue expression 
Ag {x^^w/b,x^'^w/b,x^b/w,x^'^b/w;x'^^) 



{x^w/b, x^'^w/b, x'^^b/w, x^'^b/w); x^ 



to obtain (with b — —1) an expression in elliptic functions of nome x^'^'^ 
Ag , E(-x^ /w)E(-x^^ /w)E(~x^^ /w)E{~x'^^w)E(-x'^^w)E(~x^'^w) 

— in 

Ao E{~x^w)E{~x^'^w)E{-x^^w)E{-x'^'^/w)E{-x^^ /w)E{-x^^/w) ' 

(A6) 



4 Mass 777-7 
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We begin with Wj — aj for j — 1, . . . , iV — 6 and wn-5 = bix^ , WN-i — b2X 
wn-3 = 632;"^^, u>jv_2 — b4X^'^, wn-1 = b5X~^^, wn = bgx^'^. Once again the 
bi{= b) are all equal. The Bcthe equation for the other roots is 



E{x^'/a) 



Eix^^a) 



N 



^ , , 6.3/5 E{xH/a)EHx'%/a) 
y ^-6 ) 56 E{x^a/b)E^{x^^alb) 

E'^{x^%la)E{x^%/a) E{x'^''a/aj)E{x^''ajla) 
^ E^{x^^a/b)E{x^^a/b) jj^ E{x^^aj/a)E{x^''a/aj)' 

In the X limit this gives 
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which leads to the expression in the various coefficients 



UJ 



and from the six Bethe equations involving 6, 



-{An-& 

UJ 



^6^3/5 



1. 



The recurrences to be solved for the auxiliary functions are 

^30a ^30a ^34a.^2r^ T^(^2s 



^7(a) = 



C„22a 26a 26a 3Ua 3Ua 34a 2r\ |'^-2s„^ 
V 6' 6' b' b' 6' b' )oo •'Ty-l' U') 



(^SSl rylOa loa 140 14a ^182. 2r'i '^„('^.4s,^^ ' 
V b' b' b' b' b' b' *^7V'^ ^) 

"■A,^38A,^38|^^42^.^2.)^^^(^2./„) 



/-^3Ub 34 b . 



(a;46 A 2,50 fe 2:50 ^ , 2-54 b 3,54 b 58 6 . a.2r) Q (a;4s/a) ' 
^ a' a' a' a' a' a' '°° ^ ' ' 

which have solution 



Ma) = 



('t.34o 38 a ^42 a 46 a 2r^, 
\ b' b' b' b' ^"^ 

('™10a ^14a 18a 22a 2r'\ 
V b' b' b' b' ^"^ 



(^2;30a 2;34a ^38 a 



_ T.42a 42 a 46 a 50 a 54 a 12s ^ 

J 



^10 a 1 



18 a ^66 a 70 a ^74 a 78 a 12s 



^7(1/0) 



C™34 6 „38b ™42 5 „46 b . „2r^ 
. '■•^ g'-^ a'^^ g'-^ g'*^ 
'(^58b ^62 b 66 i ^,70i.™2r^ 

^•^ a'-^ a'-^ a'*^ a'"^ 

/'„30b ^34 b 38 b ^42 b 90 b ^94 b 98 b ^102 b 12s- 

V , ^ a ' a ' a ' a ' a ' a ' a ' a ' 



00 



/ 54 b ^58 b ^62 b 66 b ^66 b ^70 A „74b ^78b I2s\ 
g'-^ a'-^ g'-^ g'*^ a'"^ g'*^ g'*^ a'"^ = 

Substitution into 

A ,„4 (''r22 2£ ^26™ ^30 H 34 u; 6 b_ „10b_ ^14 A. ^18A•^.2r•^ 

j,;-^ b'-^m'*^ M)'"^ M)'"^ u-'*^ 

3 ('xSH -rl02£ -rl4 2£ t-18 2£ ^22 A -r26 A t30A T.34A--r2r-^ 

\ 57 57 5' ^7 

X ^7(a;''«')e7(l/^''H, 
yields the result (with b = —1 and elliptic nome a;^^*) 

A7_ 4 E{-x^/w)E{-x^^/w)E{-x^^/w)E{-x^^ /w) 
Ao ~ ^ i;(-a;6u;)£^(-a;iOw)i;(-a;i4u;)£;(_.^i8y,) 

£;(-x42w;)£;(-x46w;)£;(-x50w)i;(-x54w) 



X 



E{-x'^^/w)E{-x'^^ /w)E{-x^° /w)E{-x^'^/w) ' 
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Table I: The integers appearing in (|) and (I) for L = 3. 



j a 

1 1, 11 

2 7, 13 

3 2, 10, 12 

4 6, 10, 14 

5 3, 9, 11, 13 

6 6, 8, 12, 14 

7 4, 8, 10, 12, 14 

8 5, 7, 9, 11, 13, 15 



Table II: The integers appearing in (0) for L — 4. 



1 6 

2 1, 7 

3 4, 8 

4 5, 7 

5 2, 6, 8 

6 4, 6, 8 

7 3, 5, 7, 9 
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Table III: The integers appearing in (^) for L — 6. 



3 


a 


1 ,1 


4 


2 


1,5 


3, 3 


3, 5 


4 


2, 4,6 



Table IV: String positions Uj and corresponding eigenvalue bands for the seven 
elementary mass excitations rrii of the dilute model in regime 2^ |^. The 
strings are in units of 7r/20. 

1 String positions Band 
~ ±2,10 w 

2 ±7 

3 ±6,10 

4 ±3, ±9 

5 ±6, ±8 

6 ±4, ±8, 10 w3 

7 ±5, ±7, ±9 w'' 



Table V: The four perturbations of the tricritical Ising model, and the objects 



im statistical mechanics to which they are 


related in the scaling 


limit. 


Perturbation 


Field 


Weight 


IRF model 


^BEG 


Leading magnetic 


<t>{2.2) 


3 

80 


Not integrable 


H 


Leading thermal 


0(1,2) 


1 

10 


Dilute A4, regime 2 


l/J 


Subleading magnetic 


0(2,1) 


7 

16 


Dilute A3, regime 1 


H3 


Subleading thermal 


0(1,3) 


3 
5 


ABF A4, regime III 


D 
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